NASA TM X-71907 


NASA TM X-71907 


NASA TECHNICAL 

MEMORANDUM 


(NASA-TK-X-71907) II MS TEMPEMTaRE-STRBSS 
DEPENDENCE OF BOFGN FIBEH DEFOBKATION (NASA) 
35 p HC $4.00 CSC! 11D 


G3/24 

FADILfe 

As'-' 


N76-21288 


Unci as 
21545 


TIIV\E-TEIVIPERATURE-STRESS DEPENDENCE 
OF BORON FIBER DEFORMATION 

by James A, DiCarlo 
Lewis Research Center 
Cleveland, Ohio 44135 

TECHNICAL PAPER to be presented at 
Fourth Conference on Composite Materials Testing 
and Design sponsored by the American 
Society for Testing and Materials 
Valley Forgej Pennsylvaniaj May 3-4, 1976 




E-8698 


TIME-TEMPERATURE-STRESS DEPENDENCE 
OF BORON FIBER DEFORMATION 
J. A. DiCarlo^ 

Lewis Research Center 
ABSTRACT 

Flexural stress relaxation (FSR) and flexural internal friction (FIF) 
techniques were employed to measure the time-dependent deformation of 
boron fibers from -190 to 800° C, The principal specimens were 203 jJ-m 
(8 mil) diameter fibers commercially produced by chemical vapor deposi- 
tion (CVD) on a 13 M m tungsten substrate. The observation of complete 
creep strain recovery with time and temperature indicated that CVD boron 
fibers deform flexurally as anelastic solids with no plastic component. 

The mechanical relaxation processes responsible for the anelasticity were 
found to be controlled by a continuous distribution of thermally activated 
relaxation times t exp (Q/kT), Excellent fit of the FSR data with 
FIF data was obtainable across more than seven time decades by assuming 
that all processes have the same In -33.1) but are continuously and 
broadly distributed in energy Q, The incorporation of the derived distri- 
bution parameters into anelastic theory yielded simple analytical equations 
to describe the time-temperature-stress dependence of CVD boron fiber 
deformation at low strains. In order to extend the flexural results to high 
strain and to conditions of axially applied stress ^ an analysis of fiber axial 
deformation was made using data from the literature and from our own 
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fracture measurements. For the latter study etched fibers were employed 
to ensure core-initiated fracture, A model of a ’^composite’ ^ fiber with 
an anelastic sheath and elastic core was assumed in order to convert 
fracture stress versus time-temperature data to anelastic relations de- 
scriptive of axial deformation. These results in combination with litera- 
ture data confirm the flexural result that up to at least 800^ C all nonelas- 
tic axial behavior of CVD boron fibers can be explained solely by anelastic 
mechanisms. Diameter^ core coating, etchir^, and manufacturing source 
appear to negligibly affect boron sheath anelasticity. With a simple em- 
pirical model to account for observed stress effects, final analytical equa- 
tions and curves are presented to describe total axial deformation strain. 
These relations are further developed to demonstrate the significant effects 
of anelasticity on such fiber/composite properties as modulus, creep, 
creep recovery, stress relaxation, and dampiri^ capacity. Finally, under 
the elastic- cor e/aneias tic -sheath model, modern boron fibers on tungsten 
substrates are shown to have predictable fracture stresses for time- 
temperature conditions ranging from impact to long time stress rupture. 
Possible techniques for altering these stresses are discussed, 

INTRODUCTION 

In the design and analysis of metal matrix composites reinforced by 
chemically vapor-deposited (CVD) boron fibers, the deformation character 
of the fiber is generally assumed to be elastic , Support for the elastic 
approximation arises from observations of small fiber axial creep near 
the maximum use temperatures of aluminum and titanium matrix compos- 
ites [1, 2], This creep is usually accepted as predominantly plastic in 
origin. Recently, however, flexural studies have indicated only anelastic 
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deformation in CVD boron fibers with no evidence of plasticity at com- 
posite use temperatures and above [3, 4]. 

Our objectives in this paper are twofold. First, we will demonstrate 
that up to at least 800^ C (1412^ F) all nonelastic behavior observed during 
the axial deformation of CVD boron fibers can be explained solely by 
anelastic mechanisms. This will be accomplished by first deriving anelas- 
tic parameters and equations descriptive of our low strain flexural relaxa- 
tion data. These anelastic relations combined with a simple empirical 
stress model will then be shown to describe well the ti me- temperature- 
stress dependence of axial deformation data taken from both published 
reports and our own fiber fracture experiments. Our second objective is 
to demonstrate that, although creep strains are small, boron fiber anelas- 
ticity can produce significant mechanical effects which would otherwise be 
neglected under the elastic approximation. This will be accomplished in 
the Discussion Section by considering the design implications of anelasticity 
for such fiber/composite properties as modulus, creep, creep recovery, 
stress relaxation, and damping capacity. In addition, the role of boron 
creep in the fracture of modern fibers vapor deposited on tur^sten wire 
substrates will be considered. These particular fibers are suggested to 
have controllable and predictable time-temperature- dependent failure modes, 

ANELASTICITY 

The theory describing anelastic phenomena is discussed quite well 
elsewhere [5], Our only purpose here is to introduce briefly those con- 
cepts and relations necessary for a practical understanding of boron defor- 
mation. 

Upon application of a unidirectional stress a to an anelastic solid at 
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time t = 0 and temperature one obtains a total deformation strain 

6p(t,T,a) = eg(cr,T) + eJt,T,a). (1) 

-1 

Here = crJ^(T) = ct [E^(T)]“ is the time-independent elastic strain 
which depends on T only through the unrelaxed compliance or 
Young's modulus E^o The anelastic strain is zero at t == 0, but 
develops with time^ temperature, and stress o The three deformation 
properties which characterize and distinguish it from plastic behavior 

A 

are 

L Ldnearity: €_ is directly proportional to a, 

IL Equilibrium: After passage of sufficient time, reaches or 
relaxes to a unique equilibrium value, 

III, Recoverability: Upon removal of a, the developed completely 
disappears at a rate which is time and temperature dependent. 

From property I it is important to realize that an anelastic solid will 
creep at any stress level. This behavior can be considered to give rise 
to a time- temperature-dependent compliance 

J(t, T) = T, a)/a = J^(T)[A(t, T)] (2) 

where the stress-independent anelastic strain function A is defined here as 


A(t, T) 



(3) 


It follows, therefore, that accurate values of J^(T) can only be obtained 
near times of t - 0. In addition, for a stress -relaxation experiment in 
which a constant strain e is applied, property I implies that a(t, T, c) is 
linear with e. Thus, as with the compliance, anelasticity gives rise to 
a time- temperature -dependent modulus 

E(t, T) = E^(T)[m(t, T)] 


( 4 ) 



The strain-independent modulus ratio m defined here by Eq, (4) decreases 
with time from a value of unity at t = 0. In general, A(t)m(t) 1. 

Phenomenologically anelasticity can be described as the stress-induced 
relaxation of some micro structural entity of type i from one internal posi- 
tion to another with a relaxation time Continual motion is inhibited 
by the buildup of local elastic stresses which oppose the applied stress. 

Thus the resulting anelastic strain follows a time dependence 

e^(t, T, 0) = €g(T, ff)A.[l - exp(-t/T^)] (5) 

where the relaxation strength A. is a measure of the number and extent 
of type i relaxations. Usually, because of internal variations, the relax- 
ing entities are continuously distributed in relaxation time. One can 
account for this analytically by a continuous function x(ln r) which is de- 
fined through the total relaxation strength A; i. e, , 


A = 


X' 


X(ln t)d(ln'r) 


( 6 ) 


If the T distribution is broad (as in CVD boron), the bracketed term 
in Eq, (5) can be approximated by unity for < t and zero for r. > t. 
Thus for the broad distribution approximation (hereafter BDA), Eqs, (3), 
(5), and (6) yield 


A(t, T) = 



(7) 


thereby permitting determination of total deformation from Eq. (2) 
and a knowledge of x(ln r). Temperature variation appears in Eq. (7) pri- 
marily because the relaxation processes in boron and many other systems 
are controlled by thermally activated relaxation times; i. e. , 

r = exp(Q/kT). 


( 8 ) 
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The parameter Q is the height of the energy barrier for relaxation, 

Tq = (27rf^)"^ where f^ is the number of times per second the relaxing 
entity attempts to cross the barrier, and k is Boltzmann^s constant. 
Clearly distributions in r can arise from distributions in or in Q, 
or in both parameters. 

FLEXURAL DEFORMATION 


The unusually large anelasticity of CVD boron fibers is demonstrated 
in our recent measurements of flexural stress relaxation (FSR) and 
flexural internal friction (FIF) [4]. The principal specimens were 203 }im 
(8 mil) diameter fibers commercially supplied by Avco Systems Division. 
These fibers were produced in a single stage reactor by the hydrogen 
reduction of boron trichloride on a 13 jum (0. 5 mil) tungsten wire 
substrate held near 1300^ C (2372^ F). The FSR and FIF data of 
Ref. 4 will now be reexamined for the purpose of deriving the fiber 
anelastic strain function A(t, T). 

For the FSR data shown in Fig. 1, the stress or modulus ratio was 
determined from 


ff(t, T, e ) e (t, T, e ) 

m(t, T, e) = o_ ^ ^ . 

a(0, T, 


As indicated in Fig. 1, values of the ratio were measured optically 

3. V 

by the ratio X /X for cantilevered specimens held at low strain by end 
a o 

displacement X^ and by the ratio R^/R., for specimens held at high 

strain by constant radius of curvature R . The anelastic deformations 

o 

X and R were produced by heat treating the fibers in vacuum at tem- 
a a 

perature T for a time t of 1.0 hours, quenchir^, and removing the 
applied deformations X^ and R^, respectively. The conclusion that the 



FSR data were entirely due to anelastic sources was based on the observa- 
tion that ail flexural deformations X and R could be completely re - 
moved by reheating the fibers under zero applied deformation (property of 
complete recoverability). In addition^ the low strain m data of curve a 
was independent of end displacement and diameters between 102 and 203 i^tm 
( 4 and 8 mil). As will be discussed, the increased stress relaxation at 
low temperatures for curve b was presumably caused by a stress -induced 
activation of anelastic processes which were inhibited from operating at 
low strains by some unknown micros true tural locking mechanism. 

The extended temperature range of the FSR indicated that the deforma- 
tion strains were produced by thermally activated relaxation processes 
broadly distributed in controlling parameters. To determine whether this 
distribution was in Q, or in both, the low strain FSR data were sup- 
plemented with low strain (~10"^) FIF data from -190 to 800^ C (-310 to 
1472 ® F), The log decrement 6 was measured during free decay of 
cantilevered fibers vibrating at their fundamental and first overtone reso- 
nant modes [6]. Operating at resonance allows one also to determine the 
temperature and frequency dependence of the fiber’s Young’s modulus. 

At 22® C (72® F) and frequencies from 10^ to 10^ Hz, all modulus data 
fell within 393 ±14 GN/m^ (57 ±2 Msi). Figure 2 presents typical high 
temperature FIF data at two modes of the same specimen. The best fit 
by far to both the FIF and FSR data was obtained by assuming all respon- 
sible relaxation processes had the same but were continuously and 
broadly distributed in Q. In this situation it is convenient to introduce 
the enery distribution B(q)dq which is defined through the total relaxa- 
tion strength 



( 10 ) 




B(q)dq = / x(ln r)d(ln r) 
^InTo 


Here q = Q/k^ and the lower limit on the x integral has been 
replaced by the In of the shortest relaxation time Employing the 
broad distribution approximation (BDA) plus EqSo (7)^ (8), and (10), one 
derives the anelastic strain function 


A(t,T)-A(q) = 




B(q)dq 


( 11 ) 


where now the energy parameter 

q = (In.t - ln.TQ)T(10'^) (12) 

contains both the time (sec) and temperature (^ K) dependence. 

The best fit situation was determined by adjusting in in the q 
parameter until the A’s calculated from all sets of low strain flexural 
data agreed. For the FIF data BDA allows one to convert frequency f to 
time t = (27j-f)“^ and to calculate A from 


r 


aFIF/ X 1 
A (q) = exp} 




(13) 


.-1 


I n^ 2(10^)6{t,T) J 

Mo 

L ^ ^ Jt=(2rf)’ 

For In -33c 1 all FIF data yielded the A(q) solid line of Fig. 3. Be- 
cause internal friction was still finite below -190*^ C (-310^ F), an 
extrapolation of the 5 data to 0^ K was necessary for use of Eq. 13. 
Although the relaxation processes responsible for 6 below q = 10 make 
a negligible contribution to A , their existence implies that some 
anelastic deformation occurs in CVD boron fibers even at room tempera- 
ture. Thus the assumption inherent in the FSR measurements that 

Si 

after quenching was stable at 22® C (72® F) was not strictly correct. To 
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account for this effect in calculating A for the low strain FSR data, the 

BDA result Am = 1 was employed so that 

AFSR(q) ^ T)/m(t, T) = A(q^)/m(q) (14) 

Here A(q^) is a measure of the anelastic strain lost between removal 

at t = 0 and the FSR measurement at t^ and T = 295^ K, For q^^ ^ 12, 
BIT*’ 

^ (qj^) = lo025, and the low strain m data of Figo 1, one obtains the 

open circle FSR results of Fig, 3=. In the q determination, in was 

taken as -33 o 1 and In t was replaced by In t+ 0, 6 as discussed in Ref, 7 

to account for both the broadness and quasistatic nature of the FSR, 

The excellent agreement between A and A below q = 25 

clearly supports the model in which relaxation processes are broadly 

distributed in Q but have essentially the same In The A agreement 

is especially dramatic when one considers that some of the data differ in 

time by almost eight orders of magnitude. The deviation in the A curves 

above q = 25 has previously been interpreted [4] as evidence of another 

set of relaxation processes (Type Ib) still distributed in Q but with 

In -26o However, as will be developed shortly, it appears that a 

more appropriate explanation for the A difference at high q is the dif- 

-3 

ference in applied strain between the FSR data (e = 10" ) and the FIF 
data (e ^ 10"^). 

Althoiigh one could now determine the low stress/low q total flex- 
ural deformation of CVD boron fibers by Eq, (2) and the A(q) of Fig, 3, 
it was clear from curve b of Fig. 1 that a significant charge in A(q) 
occurred at high stress levels. Because investigation of this stress 
effect was complicated by the existence of stress gradients in flexed 
fibers, we undertook a study of the axial deformation of CVD boron fibers 



both by our own fracture experiments and by an anelastic analysis of axial 
data in the literature* It was our hope to understand not only the stress 
dependence (at least empirically) but also the extent to which anelastic 
processes contribute to total axial deformation* Included in the axial 
data, for example, would be the relaxation effects of the core region which 
was subjected to negligible stress during the flexural measurements, 

AXIAL DEFORMATION 

Approach , - A new indirect approach was developed to determine 
the time-temperature- stress dependence of boron fiber axial deforma- 
tion, This approach is based on the recent tensile results of Smith [8] 
who measured at 22^ C (72^ F) the ultimate tensile strei^th (UTS) of 
203 lira {8 mil) fibers from the same commercial run that provided the 
previous flexural specimens. Smith found that after etching a few (im 
off the surface of as -received fibers, essentially all cases of fracture 
could be explained by crack initiation within the region of the tungsten 
boride core. Apparently after surface removal, those flaws that remained 
in the boron sheath could not initiate fracture at the level of the core flaw 
strength. Smith observed that the UTS had a coefficient of variability of 
less than five percent and increased monotonically as etched fiber diam- 

9 

eter decreased, reaching values as high as 5.7 GN/m (820 ksi) near 
102 m (4 mil). 

The above UTS results of Smith suggest that the etched fiber frac- 
tures whenever its total axial deformation €p(t, T, a) becomes equal to the 

Q 

core fracture strain Assuming the core region behaves as a brittle 
elastic material (cf, appendix I), our approach to axial deformation was 
simply to measure the fiber fracture stress cTp(t, T) required for the 
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c 

boron sheath to attain Neglecting the stress contribution of the 
small 17 jLtm (0. 67 mil) core and assuming only q-type anelastic behavior, 


one can then calculate the axial strain function A(q) from Eq« (2): 


A(q) = 


E^(T)4(d) r E^(T) ap(o, d) 


(7p(q, d) E^(0° K) ap(q, d) 


where d is the diameter of the etched fiber. Smithes observation of a 


low coefficient of variability in suggested that sufficiently accurate 
d) data could be obtained. Possible sources of temperature varia- 
tion in ep(d) or ap(o, d) are discussed in appendix I. Finally, because 
modulus is directly proportional to the square of the resonant frequency f, 
the unrelaxed modulus ratio in Eq, (15) was determined from f(T) as 
measured in FIF (cf. Fig, 2); i, e, , 

_ p(T)1^ T) 

Ey(o) [f(o)J 1:+X(T) 

The slight effect of thermal expansion on f(T) is contained in the param- 
eter X(T) (cf, appendix II), The calculated E^(T) is plotted in Fig, 4 in 
terms of the ratio E^(T)/E^ (295^ K), Etching was found to have no 
effect on E^(T). 

Procedure and results , - Etched specimens for the fracture measure- 
ments were prepared in the following manner. As -received 203 /im (8 mil) 
fibers of length {21 ~ 2, 5) cm were inserted vertically into a I (>10) cm 
deep etchant solution of two parts nitric acid and one part water at 100^ C. 
Upon removal of 30 jum (1, 2 mils) off the diameter in about five minutes, 
the fiber was reversed in the solution until 30 4 m was etched off the diam- 
eter of the other end. After a cleaning in distilled water followed by 
methanol, the final dimensions of the etched fiber consisted of a 2. 5 cm 



(1,0 in, ) long center section with 143 /J-m (5,6 mil) diameter and (I -2,5) 
cm long end sections with 173 jU-m (6, 8 mil) diameter. This geometry- 
insured fracture data from the 2,5 cm (1,0 in,) center gage section. 

For calculation of A(q), it was necessary to determine o-jp(Oj d)^ the 
fracture stress at d = 143 jitm and q or T = 0, Because test tempera- 
tures at absolute zero were not possible, was measured near -195 ^ C 
(-319^ F) by surrounding the gage lei^th of the etched fibers with liquid 
nitrogen and pulling the specimen to fracture in a standard tensile test. 

For such a test, the BDA allows one to calculate q simply from the time 
of stress application. For q = 3,0, the average UTS ± standard deviation 
of 10 specimens was 5, 54 ±0,30 GN/m^ (804 ±43 ksi). When this result 
is compared to Smith- s q = 11,8 values of 5, 03 ±0, 16 GN/m^ (730 ±23 ksi), 
the significant effect of boron sheath creep as low as room temperature 
becomes obvious. Assuming at low q a direct proportionality between 
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A and q, one calculates an extrapolated average ^^(o, d) of 5, 75 GN/m 
(834 ksi). As discussed in appendix I, this Op value should increase 
slightly as test temperatures exceed the etch temperature of 100^ C (212^ F). 

For the high temperature fracture tests, the etched fibers were in- 
serted into a vertical furnace heated by a 7, 6 cm (3 in. ) long tantalum 
element with a 1,3 cm (0, 5 in, ) bore. After the fiber ends were mechani- 
cal clamped with the gage section centered in the furnace, the lower 
fiber end was loaded to a stress level by the addition of a known 

weight. The furnace chamber including the load was closed and evacuated 
“ 6 

to less than 3x10 mm Hg. Heat was then applied to provide a specimen 
warmup rate of about 1, 0° C per sec. Temperature was measured with 
a chromel-alumel thermocouple located '^2 mm from specimen center. 



Thermal gradients were less than 10^ C along the gage length. Furnace 
shutdown was automatic upon fiber fracture. Diameters of the fractured 
ends were measured optically in two perpendicular directions with an 
accuracy of 0.5 fxm (0.02 mils). 

Fracture stress data for a few fibers etched to a 143 Mm diameter 

are shown in Fig, 5 as a function of fracture temperature T^. Optical 

microscopy indicated that in all cases fracture was initiated by flaws 

within the core region. To determine q for each high temperature data 

point, Eq. (12) was maximized under the condition that both t and T 

were varying during warmup rate y. It can be shown that as long as 
0 

1 < Tp/y < 10^ sec, the proper time and temperature for the q param- 
eter are 0.01 Tj^, and T^, respectively. Thus the closed circle points 
of Fig. 5 can be considered as stress -rupture data where the failure 
life is about 10 sec. 

With consideration of possible temperature variation in o-p,(0, d) as 
discussed in appendix I, the high temperature fracture stress data with an 
assumed five percent coefficient of variability were inserted into Eq. (15) 
to determine the axial anelastic strain function A(q). The results are 
shown as curve Ajj in Fig. 6. It is assumed a priori that: the relaxation 
distribution parameters determined for q < 2& apply for all data indi- 
cated. Also included in this plot is the low strain Aj^ function from 
Fig. 3, Rather than distinguishing stress modes, the subscripts H (high) 
and L (low) indicate the qualitative stress level at which each A data set 
was taken. 

To evaluate the validity of our approach, the anelastic A(q) equations 
were applied to axial deformation data within the literature. Because 


anelastic creep occurs at all stress levels, accurate A derivations re- 
quire total deformation versus time curves rather than approximate creep 
strain curves. This requirement is met in the 102 jiim (4 mil) fiber data 
of Rose and Stokes at 538^ C (1000^ F) [9], When their time-dependent 
results at three stress levels are converted to A(q) by Eq, (2), the short 
solid lines near q = 35 in Fig, 6 are obtained. Regarding measurements 
similar to our own fracture tests, Veltri and Galasso [10] have determined 
the UTS ratio (jp(T)/crp (295® K) for 102 p-m (4 mil) boron and silicon 
carbide- coated boron fibers. Assuming a temperature -independent fracture 
strain, a few second tensile test, and A (295® K) = 1. 15, one can average 
their UTS ratios at each test temperature (except where surface contamina- 
tion is obvious) and convert the data to A(q) by Eq, (15), The calculated 
points and their variation are shown by the cross data in Fig, 6, It is seen 

that both the high stress creep and UTS axial data yield A(q) values in good 
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agreement with Aj^ despite a variation in stress from 1. 4 to 5, 5 GN/m 
(200 to 800 ksi). This stress independence of Aj^ indicates a linear rela- 
tionship between strain and stress, thereby implying the validity of our 
original assumption of complete anelasticity in boron fiber axial deforma- 
tion, In addition, the equivalence of the various sets of data suggests that 
above ~1, 7 GN/m (250 ksi) the axial Ajj is independent of diameter, 
manufacturing source, fiber coating, and etching. 

Regarding axial deformation at low stress, Ericks en [11] recently 
found that the empirical relation 

e£)(t) = C(o-) + Q?(a)log t 

adequately described the total anelastic deformation of 102 /xm (4 mil) 
silicon carbide-coated boron fibers at 27® C (81® F) and t(h) from 0, 01 



100 « From the anelastic relations in Eqs, (3), (11), and (12), it follows 
that at constant T 

a(q, ff) = 2.3(10'^)egTB(q, cr) 

Thus Ericksen's ot parameter can be employed to determine the stress 
dependence of B(q), the number of relaxation processes at q. For 
example, taking q^ = 12,6 for one hour creep at 300^ K, one calculates 
the B(q^) versus a plot of Fig. 7, The closed circles are from 
Ericksen’s fiber data whereas the squares are derived from his B/Al 
composite data assuming no matrix creep. Also included in Fig. 7 are 
the B(q^) from our own low strain FIF and FSR data (open circles). Com- 
parison of Fig. 7 with the low q behavior of Fig. 6 indicates that the 

B(q) responsible for A^^ is small at low stress. However, as stress 

2 

exceeds some critical ’’breakaway” stress a* of about 1,0 GN/m 
(150 ksi), B begins to rise rapidly, eventually saturating to produce the 
stress -independent Ajj, The maximum A(q^) corresponding to Ericksen’s 
data is plotted in Fig, 6 under the low q assumption of a direct propor- 
tionality between A and q (i, e. , B independent of q). 

It is of interest to apply this concept of a breakaway stress a* to 
the Fig, 1 FSR data. For example, since the surface stress for curve b 
was greater than 2.2 GN/m (320 ksi), the maximum number of relaxa- 
tion processes should have been operating so that at low q the function 
Ajj is appropriate. Applying the Eq, (14) correction Ajj(q^) = 1. 15 to 
the high strain FSR data, one obtains the Fig. 6 curve labeled Ag (open 
squares). Although agreement at low q was forced by the cor- 

rection, the equivalence above q = 30 between Ag and ^B is a signifi- 
cant result with many fundamental implications. First, it is another 
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result suggesting that boron fiber axial deformation is completely anelas- 
tiCo Second, the core region which was not stressed during flexure 
appears to contribute negligibly to axial deformation<, Third, as also 
indicated in Fig, 7, fiber deformation is essentially independent of the 
mode of stress application. Furthermore, if the same correc- 

tion is applied to the low strain FSR data above q = 30, the calculated A 
points in Fig, 6 (open circles) are also in agreement with and Ag, 
This result suggests that the stress dependence of boron fiber anelasticity 
disappears near q = 35. In terms of the breakaway stress model, one 
would conclude that a* tends to zero above q = 35, The Fig, 3 results 
further suggest that a* is large up to q = 20 but then starts to diminish, 
resulting in the strain- dependent deviation near q = 25. One can incor- 
porate these concepts into the approximate q dependence for a* shown 
schematically in the Fig, 6 insert. 

On the practical side, a stress model which incorporates an abrupt 
transition from Aj^ to at a certain stress level or* is certainly a 
crude approximation, especially when one considers the residual stress 
pattern in CVD boron fibers [12], Nevertheless, the Fig, 7 results do 
indicate that such a model can serve as a first approximation for deciding 
which A function is operating under certain stress conditions. Ericks en's 
data and our own FSR data suggest a cr* value somewhere near 

XjlIcLX 
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1,0 GN/m (150 ksi) for q < 20, As temperature or q is raised, a* 
should diminish somewhat in the manner indicated in the Fig, 6 insert. 

On the basic side, the above stress behavior suggests a microstruc- 
ture model in which the application of high stress or high temperature 
loosens low q relaxation processes from internal locks so that they can 



operate and be observed, Firle in his internal friction study of 102 jum 
(4 mil) fibers [13 ] came to a similar conclusion by suggesting a model of 
stress-induced dislocation breakaway. More recent work, however, has 
found no evidence of dislocations in CVD boron fibers [14]. Therefore, 
in keeping with our model previously developed [4], we suggest at the 
present time that the relaxing entities are boron icosohedron and the in- 
ternal locks are unknown atomic imperfections quenched in after deposi- 
tion. The dropoff of a* with q would then indicate some form of ther- 
mally assisted breakaway, perhaps caused by thermally induced mobility 
of the internal locks. 

DISCUSSION 

Although many basic microstructural implications can be drawn from 
the previous results [4], our only purpose here is to discuss briefly some 
important practical effects of an elasticity on the mechanical properties of 
CVD boron fibers and their composites. 

Summarizing the results of the anelastic approach, one can com- 
pletely describe the time-temperature-stress dependence of boron fiber 
axial deformation by the simple relation 

ej5(t, T, a) = A (q) (17) 

where 

q = (In t + C)T(10‘^) 

and 

\ = Aj^(q), a < a*(q) 

\ = Ajj(q), ff > ff*(q). 

The Aj^, Ajj, and cr* are shown as functions of q in Fig. 6. The value 
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of a* „ is estimated to be about 1,0 GN/m (150 ksi). For quasistatic 
max 

tests such as tensile, creep, and stress -relaxation, t (sec) is the time of 
stress (strain) application and C = 33,7o For dynamic measurements at 
frequency f, a is the stress amplitude, t = (2?7f)~^, and C = 33,1. E^(T) 
can be derived from Fig, 4 with a (295^ K) value of 390 GN/m^ (57 Msi), 
For flexural tests, Ajj is replaced by Ag of Fig, 6, and and a can 
be taken as the surface strain and stress, respectively. 

With this analytical base, let us now consider the anelastic effects on 
some important fiber properties. 

Modulus , - The anelastic fact that the total deformation strain is di- 
rectly proportional to a leads to the significant practical result that 

E(t, T) = aAp(q) = EJT)/A„(q) (18) ■ 

Thus for Young ^s modulus measurement, the length of time during stress 
application can be critical. For example, in the low stress regime near 
room temperature, A^ = A^^ = 1 so that one can effectively measure E^ 
from a stress -strain curve. However, at hi^h stress or high temperature, 
A^(q) > 1 so that stress -strain slopes gradually decrease durii^ the time 
of test, thereby yielding an average modulus less than E^(T). Reported 
data which shows the E for boron decreasing rapidly with temperature 
and time [15, 16] are clear evidence of the anelastic effect. 

This modulus behavior may become quite significant for fibers sub- 
jected to dynamic stress conditions at high temperature. For example, 

for low stress amplitude at 482^ C (900^ F) and f = 10^ Hz, published 

2 

tensile test modulus data [15] yields E = 240 GN/m (35 Msi), whereas 
the dynamic modulus using the above analysis is 350 GN/m (51 Msi). 
Increasii^ the frequency, of course, lowers the stress application time 
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and raises the effective modulus toward E^(T) = 370 GN/m"^ (54 Msi), 

Creep . - Because anelastic solids creep at all stress levels, it is 
difficult in a fiber creep test to separate the elastic strain from the time- 
dependent anelastic strain without quenching to lower temperature and 
removing the stress. Therefore, in discussing creep it is more appro- 
priate to consider the total deformation as given by Eq. (17). Because 
of the In t dependence in the q parameter, at a given temperature 
will increase rapidly at short time and then level off somewhat as time 
progresses. Although it may be convenient to label this behavior as 
primary and secondary creep, due to a difference in mechanisms one 
should not equate it with similar creep behavior observed in metals. 

Regardir^ creep rate, differentiation of Eq. (17) yields 


T, ff) = 

^ Eu(T)t [dqj 


Thus in certain q regions where A is fairly constant (cf. Fig. 6), the 
"secondary'^ creep rate will vary inversely with time, until q values 
are reached where A increases again. Fundamentally it should be pointed 
out that macroscopic equilibrium or exhaustion of the anelastic processes 
is not observed in boron fibers because these processes are continuously 
distributed in energy or q. That is, when one process reaches micro- 
scopic equilibrium, another one at higher energy is activated so that creep 
strain saturation is never reached (at least to q = 45). 


Finally, in those situations where different stresses a. are applied 
to a boron fiber at different times t^ the Boltzmann superposition prin- 
ciple allows one to determine e^(t) by summing the deformations from 
each stress as if each were acting alone [5]; i. e. , 
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e^it, T) = [E^(T)]-l^ffjA.(t. - t, T) (20) 

For continuous stress application^ the summation is replaced by an inte- 
gral^ as, for example, during a tensile test or during composite creep 
where load transfer to the fibers occurs during matrix stress relaxation. 
For these situations, because of the In t dependence one introduces 
negligible error by employing Eq. (17) and assuming t is the time from 
zero to maximum stress a. 

Creep rec ove ry. ~ When the stress is removed from a solid which 
has crept anelastically, the creep strain gradually disappears with time 
and temperature. This behavior called the elastic aftereffect or creep 
recovery can be analytically described by use of the superposition prin- 
ciple, For example, if creep occurred durir^ application of stress 
for time t^ at temperature T^, then when is removed at t = 0 and 
T, the creep strain recovers according to the relation 

T, a^) = a^[EJT)]-l[A^(q*) - A^(q)] (21) 

where 

q* = qp + (lO'^)T ln[l + (t/t*)] 

and 

In t* = 10®(qyT) - C 

In this derivation the slight temperature dependence of E^ has been 
neglected. 

The property of creep recovery may at times be a significant consid- 
eration in the design and use of metal matrix composites. For example, 
for those composites formed by stressing boron fibers at high tempera- 
ture, fiber recovery w'ill be small as long as the inequality A(q) < A(q*) is 
maintained. However, any high temperature excursions during use condi- 



tions could induce dimensional changes which would be unexpected under 
the assumption of a plastic fiber » Alternatively, recoverability might be 
employed to remove by proper heat treatment any undesired composite 
creep developed during use^ Obviously this technique requires matrices 
with near zero creep resistance at the treatment temperature. Along 
these lines it should be understood that to date complete strain recovera- 
bility has only been accurately studied under conditions of high vacuum. 
The possibility may exist that high temperature air and matrix environ- 
ments produce enov^h contamination to lock some of the micros tructural 
mechanisms in the relaxed state, thereby preventing total recovery. 

Stress relaxation , - Because relaxation processes are broadly dis- 
tributed in energy, boron fiber stress relaxation as a function of time, 
temperature, and strain is directly determinable from the anelastic strain 
function; i. e„ , 


a(t, T, e) ::: E^(T)6/A^(q). (22) 

Damping capacity , - The damping capacity ip of 3 l vibrating system 
is the ratio of the energy loss per cycle to the total stored energy. Since 
the log decrement 5 = i///2, it follows from Eq, (13) that for a boron 
fiber 


Mf, T) = 


7T 1 


10 ' 


dA 


n 

a 


A^dq 


(23) 


.^-l 


t = (27rf)‘ 

Figure 2, which plots the low stress 6(T) near 1 kHz indicates a rapid 
rise in i//(T) with significantly large values = 24 percent) within the 

use- temperature range of metal matrix composites. Although high stress 
6(T) data has yet to be taken, the low q slope of Ajj in comparison to 
Aj^ siiggests that high vibrational stresses should induce even larger \p 
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RTOUCIBlLirY OF 


ORIGIN B.I' -Gii: rui/jri. 

at lower temperatures (cf. Peak Ia**[4]). Obviously this behavior should 
be considered in vibrating composite systems, especially since the contri- 
bution from the fiber or matrix to the composite is directly proportional 
to the component's stored energy or modulus. For example, neglecting 
the interface and other sources of energy loss, one can write for axial vi- 


bration that 


comp 


(VE)fiber + 


where V is the volume fraction. Thus for a 50 vol percent boron/ 
aluminum composite, over 80 percent of is contributed by 

Clearly if boron fibers were assumed elastic with no dampii^ capacity, 
a significant mechanical property would be neglected. 

Fracture . - The temperature- dependent fracture results of this study 
coupled with Smith’s UTS results at 22^ C (72° F) [8] confirm the impor- 
tant role of the core flaw in the fracture behavior of modern boron fibers 


vapor- deposited on a tungsten wire substrate. The model that emerges is 
that of a ’’composite'^ fiber consisting of an elastic tungsten-boride core 
within an anelastic boron sheath. When etching removes surface flaws 
from commercial 203 fi m (S mil) as-received fibers, essentially no flaws 
remain in the sheath of strength less than that of the core flaw. Thus frac- 
ture in these fibers is produced by the elastic extension of the core region 
to its fracture strain as modified by internal residual stresses (cf. appen- 
dix I). The situation is not as clear in 102 jLtm (4 mil) and 143 jum (5. 6 mil) 
fibers where Smith has observed a larger variation in UTS after etching. 
This might be explained either by flaws in the boron sheath or more likely 
by incomplete core development due to reduced deposition times. 
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Fracture of etched fibers occurs ^ therefore, at the nearly constant 
fracture strain of the tungsten boride corec Fracture stress, on the other 
hand, decreases with time and temperature in a manner dependent on the 


anelastic axial creep of the boron sheatho Its value can be simply derived 


from Eq. (17) as that stress required to deform the fiber to the core frac- 


ture strain e^; io e» 


Op(t, T, d) = 

^ Ajj(q) 


E„(T)e^(d) r E ,(T)1 rAjj(qQ) 


E^(To) Ajj(q) 




Here d is the diameter of the etched fiber and stress contribution of the 


17 jum core has been neglected„ This result suggests that given a fiber 
with a core-initiated fracture stress d) at a certain value, 

one can now predict from its fracture stress at any other q, provided, 
of course, that the particular use conditions do not alter or initiate 
flaws or lower strength than the core. The practical significance is quite 
dramatic when one considers that all types of time- dependent stress appli- 
cations are included, ranging from impact tests through tensile tests to 
long time stress -rupture tests. For example, if the average UTS of a 
group of fibers is 2. 76 GN/m (400 ksi) at room temperature (Ajj = 1. 15), 

their average fracture stress at 300^ C (572° F) will be as high as 
2 

2.4 GN/m (350 ksi) for impact tests of 1 msec duration (Ajj — 1»3), but 
will decrease to 1.7 GN/m (240 ksi) for 1000 hour stress-rupture tests 
(Ajj = 1.9). Clearly the Eq. (25) calculation can be applied directly to 
unidirectional composites and to axial fracture conditions in which any two 
of the three variables time, temperature, and stress are known. 

There are two other important implications of the core fracture model 
that should be considered. On the negative side, boron fibers allowed to 
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creep at high temperature under stress and then lowered in temperature 
should have a reduced UTS at room temperature. This behavior is ex- 
plained by the fact that a certain tensile creep strain is quenched into the 
sheath and core, thereby requiring less elastic strain to achieve core 
fracture at the lower temperatures. On the positive side, however, it may 
be possible to creep boron fibers in compression at high temperature, as, 
for example, durii^ composite production. This treatment should quench 
in a compressive strain in the core, resulting in a higher UTS at room 
temperature. One calculates that if a fiber were subjected to a compres- 
sive axial stress a for one hour at 500^ C (932^ F) or one second at 
688^ C (1270^ F), its room temperature UTS would increase by ~al 

CONCLUSIONS 

Due to the analytical nature of an elastic deformation, it is now pos- 
sible to predict more accurately the mechanical behavior of CVD boron 
fibers both as free fibers as reinforcement in a composite. At low tem- 
perature and short times these fibers can be considered as elastic ma- 
terials; but as temperature or time is increased, anelastic behavior 
begins to manifest itself significantly in many nondestructive and destruc- 
tive mechanical properties. Our desire in this paper was not to discuss 
in detail the anelastic implications for a wide rai^e of design situations. 
Rather, it was our belief that a brief insight into possible effects would 
aid the composite designer in obtaining enough understanding to more 
fully exploit boron fiber anelasticity for his own composite problem. 
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APPENDIX I - CORE FRACTURE STRAIN 
For an as-received boron fiber of diameter d^ (= 203 fim) etched to 
diameter d, the fracture strain €p(d) of its tungsten boride core can be 
considered to consist of four components: 

e^i&) = ep(d^) - q) 

= [e^(0) - e^(d^) - e^(T)] - e^(d, q) (II) 

where 

ep(0) fracture strain of a core free of the boron sheath (assumed to 
be temperature independent), 

residual deposition strain produced near 1300° C by bo riding of 
tungsten substrate and by vapor deposition of boron, 
e (T) residual thermal strain due to difference in thermal expansion 

X 

between boron sheath and tungsten boride core, 
e„(d, q) residual etching strain produced by removing outer compressive 
layers of the boron sheath [e„(d^, q) = 0] 

For the particular fibers of this study, Behrendt [12] has measured near 
room temperature the sum + e^) and the etch strain (143 fim) and 
found both to be compressive (negative) . Thus the fracture strain for the 
core in the etched fiber is greater than that of a free core= 

Regarding the effects of temperature on €^(d), the constancy of 
Smithes maximum UTS data at 22° C [8] suggests that €^(d^) is not affected 
by heat treating the as-received fibers to 900° C for 1, 5 hours. (The UTS 
falloff from 600 to 900° C was recovered completely by etching, thereby 
indicating the influence of thermally- induced surface defects above 600° C. ) 
The insensitivity of ep(d^) to heat treatment further suggests that up to 
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at least 900^ C there is negligible change in and the core deforms 

elastically (io e. , no core creep). 

For e^(T), the observation of similar thermal expansion coefficients 
near room temperature [17] coupled with the temperature dependence of 
their specific heats [18] suggest throv^h Gruneisen^s equation [19] that the 
core and sheath expand similarly up to "^500^ C where the ^26^ expan- 
sion rate begins to decrease slightly, In terms of 6^(T), this behavior 
can be interpreted as a constant compressive e of about 10 until 
500^ C where begins to decrease, reaching zero at the deposition 
temperature. 

The contraction strain e^Cd, q) is the result of elastic plus anelastic 
recovery of the inner tensile regions in the boron sheath [12], Its value 
can be calculated from the superposition principle [5]: 

= o^{d)A{q)/E^(T) (12) 

where a (d) is a fictitious compressive stress produced by the removal 
of outer compressive layers , Those fibers tested below T = 373^ K in 
this study were assumed to have an value determined by d = 143 iim 
and A(q^) where q^ corresponds to the etch conditions of '^lO min, at 
100^ C, Fracture of these fibers occurs, therefore, when of the 
sheath equals €jp(d); or, from Eq, (II), when 

CTp(q, d)A(q) cr„(0, d ) ff (d)A(q') 

■J- = -E 9 5 •. (13) 

E^(T) E^(0) EJT) 

Assuming from our low q data that A(q) = 1. 15 for Smith’s UTS versus 

d results [8], we estimate 0-^(0, d^) = 5,17 GN/m^ (750 ksi) and 

2 

a (143 jum) = -0,55 GN/m (-80 ksi). For the high temperature fracture 
tests, these stress values are inserted into Eq. (13)' with appropriate 
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£ (T) correction and A(q-) replaced by A(q) to account for the additional 
contraction produced by for T > 373° K, The resulting equation is 
then employed with the experimental CTp(q, d) to determine A(q), Although 
the corrections for the thermal dependence of c and €„ are somewhat 

X jLi 

compensating and affect A only slightly, the actual occurrence of the 
effect appears to be supported by Smithes observation of higher room 
temperature UTS for etched fibers heat treated near 500° C, Obviously 
the correction for this effect could be eliminated by prior heat treatment 
above the fracture test temperatures. 



28 


APPENDIX II - BORON FIBER THERMAL EXPANSION 
Employing published data on thermal expansion over large tempera- 
ture intervals, we have calculated the boron fiber thermal expansion co- 
efficient a(T) through the Gruneisen equation [19]: 

1/2 


c ^ 

1/2 

~cl 

V 


V 

Za 


30? 


- bH 


H=0 


Here is the specific heat at constant volume (cal/mole/^ C), H is the 
internal energy (cal/mole), and b is a constant equal to ~0.01 (mole/ 
cal)^^^. Because the Debye temperature 6 of amorphous boron is ~1200^ K, 
the approximation is fairly accurate for T < 0. With the amor- 

phous Gp(T) and H(T) results of Wise et al. [20], we have adjusted o?(T) 
in the Gruneisen equation to fit the following boron fiber data: 


.-6 


a = 4.9x10'"'' per C, 300-600^ K [21] 

a = 5.4x10'® per ° C, 300-600° K [17] 

a = 8.5x10'® per ° C, 300-1300° K [17] 

a = 8.3x10'® per ° C, 300-1300° K [22] 

pT 

Our results for a and the parameter X - j adT are shown as a 

^0 

function of T in Fig. 8. 
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Figure 1, - One hour flexural stress relaxation IFSR) of 203 pm 18 mill boron fibers 
as a function of temperature and average surface strain e^. 
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Figure 2. - Flexural internal friction (FI F) and vibration frequency at two resonant 
modes of a cantilevered boron fiber. 
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Figure 6 . - Anelastic strain function of boron fibers as derived from various axial and 
flexural deformation data. Insert empirical estimate of the q dependence of the 
"breakaway" stress cf . 
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Figure 7. - Stress dependence of B(qQ>, the relaxation strength or number of an- 
elastic processes with energy q^. 
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Figure 8. - Thermal expansion parameters of boron fibers as derived from 
Gruneisen's equation and expansion data. 
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